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Методика дискретного управления электродвигателем постоянного тока* 
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Рассмотрена методика дискретного управления непрерывным объектом на примере двигателя постоянного 
тока. Предлаг. аемая методика позволяет упростить процесс выбора желаемых свойств синтезируемой СИСТЕ- 
мы управления, свести нахождение управления к ст. андартной задаче решения системы линейных уравнений 
[@ учётом ОГ! раничений-неравенств на каждом шаге дискретизации по времени. Желаемые свойства синтези- 
руемой системы задаются с помощью эталонного переходного процесса. Используется квадратичный крите- 
рий качества, харак теризующий отклонение переходного процесса от эталонного в равноотстоящие моменты 
времени. Такой подход позволяет в реальном масштабе времени управлять объектом с учётом ОГ! раничений 
как на управляющее ВОЗДЕЙСТ! вие, так и на фазовые координаты системы управления. При ЭТОМ МОЖНО ДО- 
биться отсутствия перерегулирования и колебательности. 
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Введение. Электродвигатели постоянного тока являются непрерывными объектами. В то же 
время для реализации управления всё большее применение находят средства, использующие 
представление информации в дискретной форме. Причём управление должно строиться с учётом 
нелинейных свойств объектов управления и реальных ограничений на фазовые координаты и 
управляющее воздействие. 

Для решения задач оптимального управления большое распространение получили мето- 
ды, опирающиеся на использовании минимизируемых функционалов. Существующие методы син- 
теза можно разделить на две группы: методы, использующие косвенные критерии качества, и 
методы, использующие прямые критерии качества. 

К методам синтеза оптимальных систем управления с использованием косвенных критери- 
ев качества относятся принцип максимума Понтрягина [1] и принцип оптимальности Беллмана 
[2]. Эти принципы являются мощными математически обоснованными способами решения опти- 
мизационных задач в отношении минимизации выбранного критерия качества. Синтез систем 
управления с помощью принципа максимума Понтрягина сводится к решению двухточечной крае- 
вой задачи для дифференциальных уравнений [1]. Получение аналитического выражения для 
оптимального управления в замкнутой форме связано с большими трудностями и представляет 
собой самостоятельную задачу для каждого класса объектов [3]. Трудности решения двухточеч- 
ных краевых задач стимулировали поиск разного рода прямых методов [4—6]. 

В данной статье предлагается методика сведения задачи нахождения управляющего воз- 
действия для нелинейной системы управления к решению системы линейных алгебраических 
уравнений на каждом шаге квантования по времени. При этом требуется обеспечить отсутствие 
перерегулирования, колебательности и учесть ограничения на фазовые координаты и управляю- 
щее воздействие. 

Для задания желаемых свойств переходного процесса предлагается использовать квадра- 
тичный критерий качества, характеризующий отклонение переходного процесса от эталонного в 
равноотстоящие моменты времени. 

Математическая модель объекта управления. При упрощённом исследовании электродвига- 
телей постоянного тока с использованием линеаризованных математических моделей обычно 
принимают следующие допущения: гистерезис в магнитной цепи отсутствует, то есть намагничи- 
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вающая сила однозначно определяет магнитный поток; размагничивающее действие реакции 
скомпенсировано; сопротивление и индуктивность якоря не зависят от значения тока возбужде- 
ния, тока якорной цепи и других факторов и считаются неизменными; к валу двигателя со сторо- 
ны машины приложен электромагнитный момент /М, а все моменты, обусловленные механиче- 
скими потерями, входят в момент сопротивления нагрузки М. = М. (Ё, а, ®) ; момент инерции при- 


вода не зависит от параметров его движения; управление осуществляется по одному каналу. 
Таким образом, электродвигатель постоянного тока с независимым возбуждением описы- 
вается системой дифференциальных уравнений следующего вида: 


аа 


Е ь 
ао К КК 1 

=-Ти--"®о+- М (0), 
а К УВ У 


где К — сопротивление якорной цепи; о — угол поворота вала двигателя; « — угловая ско- 
рость вала двигателя; М. (5) — момент сопротивления на валу двигателя; 7 — момент инер- 





ции, приведённый к валу двигателя; К„, К, — параметры, характеризующие электромехани- 


[о 


ческие свойства электродвигателя; и — напряжение в цепи якоря. 
Без потери общности момент нагрузки на валу двигателя будем считать равным 
М. (5) =с.®, 
где с — некоторый коэффициент. В случае нелинейного момента на валу можно осуществить 
линеаризацию на каждом шаге квантования. 
Система уравнений, описывающая электродвигатель является линейной и может быть 
представлена в матрично-векторном виде [7]: 


х(Е) = Ах(Е) + ВЦ(Е), (2) 
к 0 1 
где х(Е) = | =— вектор-столбец состояния электродвигателя; А= КК ю | — постоянная 
[© Е 
К 9 
0 
матрица параметров системы; В = К, — ПОСТОЯННЫЙ вектор-столбец параметров системы; 
К 


и(Е) — управляющее воздействие. 
Методика управления. Дискретные модели объектов строятся в предположении, что управ- 
ляющие воздействия, формируемые в ЭВМ в виде цифровых кодов, преобразуются в непрерыв- 
ные процессы с помощью экстраполяторов нулевого порядка. На практике в качестве экстраполя- 
торов нулевого порядка используются цифроаналоговые преобразователи [6]. 
Рассмотрим задачу перевода системы из начального состояния в нулевое. На управление 
и вектор состояния наложены ограничения [8]: 
И и (ЕЁ) < Их, (3) 
Хи < (Е) < Хи. (4) 
На каждом А-м интервале дискретизации по времени объект управления (2) будет описы- 
ваться уравнением [9]: 
Хх, (© =Ах,(В+Ви,, 
а решение этого уравнения в моменты дискретизации описывается линейными разностными мат- 
ричными уравнениями вида: 
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х,! =е^Тх, +Р(Гуи, , 


КТ-+Т 
где Р(Т)= | е^\7'"ЭВат; Т — период дискретизации. 


АТ 


Для решения задачи перевода вала из начального состояния х, = {© (0) = 0%; (0) = о} В 
конечное х„„ = ое = он (Е ) = ©, |} С МИниМальными потерями энергии в якоре, эталон- 


ный переходный процесс будем выбирать на основе известной процедуры построения оптималь- 
ной диаграммы тока и скорости, так чтобы переходный процесс заканчивался за 4 шага дискрети- 
зации по времени (период дискретизации Т примем равным 0,013 с). 

Определим последовательность и, таким образом, чтобы минимизировать сумму квадра- 
тов отклонений от эталонного переходного процесса в моменты квантования. Тогда критерий ка- 
чества представляется в виде: 

у5 -5‹ > тм, 
К=0 
где 5, =х, -х, — вектор ошибки размерности 2 в моменты дискретизации; х, — вектор состоя- 
ния эталонного переходного процесса размерности 2 в моменты времени (АТ. 

Тогда задача синтеза как задача о наименьших квадратах с линейными ограничениями- 
неравенствами формулируется следующим образом [10]: минимизировать 

[Е.М -=| (5) 
при условии 

СИН, (6) 
где Е, ЕЁ С, В, И— соответственно лхл7-матрица, л-вектор, /х 7 -матрица, /-вектор, п-вектор. 

Для сформулированной задачи необходимо минимизировать: 


= * т * 
УХ — Хил Хр — Хр (7) 
к=0 





с учётом ограничений (3) и (4). 
Поскольку переходный процесс заканчивается за 4 шага, выражение (7) примет вид: 


* ж 


т 
Хх, Хх, -Х, 





Х, 

з Е Хх, —Х, Хх) —Х, Е 
УЖ Ха — Хуа =: * У ти. 
К=0 Хз — Хз Хз — Хз 

Хх, -Х. Хх. -х. 


Метод наименьших квадратов для решения задачи оптимизации предполагает одновре- 
менное нахождение всех неизвестных с помощью решения одной системы линейных уравнений 


[10]. Поэтому при нахождении и, (! =0,3) необходимо в вектор неизвестных Ш включить все 
и, (1 =0, 3). 
Введём вектор неизвестных (/ = [4, ‚и, и.]”. 


Определение такого вектора неизвестных позволяет построить программное управление. 
С целью построения замкнутой системы управления после нахождения вектора неизвестных на 
первом шаге дискретизации на вход подаётся управляющее воздействие, сформированное с по- 
мощью системы базисных функций и вектора ш‚. Векторы и,,и,, и. не участвуют в формировании 


управляющего воздействия. Однако вектор и, используется для определения действия ограни- 
чений на управляющее воздействие на шаге дискретизации 4. 
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Нахождение управляющего воздействия для второго шага дискретизации по времени про- 

исходит путём решения системы уравнений с ограничениями после измерения состояния х, на 
и 

первом шаге. В этом случае вектор неизвестных равен И =| и, |. Для формирования управления 
и 

на втором шаге используется вектор и, . 

Следовательно, включение в вектор неизвестных управлений на всём периоде действия 
ограничений позволяет синтезировать оптимальную систему управления путём решения последо- 
вательности линейных алгебраических уравнений. 

Предлагаемый способ учёта ограничений позволяет гарантировать их жёсткое соблюде- 
ние в отличие от использования штрафных функций. 

Формирование матриц для решения задачи управления. Сформируем матрицу Е и вектор 
Е входящие в матричное выражение (5). Для этого найдём выражение для вектора состояния сис- 


темы и эталонного вектора состояния в моменты времени АТ +Т, где К =0,3. 

Для каждого К-го шага (к =1, 4) вычислим значения эталонного переходного процесса по 
формуле: 

х, =е“Тх.. (8) 

Переходный процесс в системе определяется выражениями: 

х.1 =е“х, +Р(Г)и,, 

х.) =е^Тх, „ +Р(Гуи, =е?^Т х, +е^"Р(Т)-и, + Р(Г)и и, 

х.з =е^х,„ +Р(Т)и,., = е3^Гх, че?" Р(Т)и, +е^ТР(Г) и, +Р(Г)и,.., 

ха =е^“х,; +Р(Гди,. =е“^" х, +е\\"Р(Т)и, +е?^"Р(Т)и,, +е^"Р(Т)-и,,, +Р(Т)и 
Согласно методу обратной задачи динамики [4, 11] подставим выражение (8) для 


(9) 











1+3 1+2 1+3 


х,.: (К =0, 3) в левую часть системы уравнений (9) вместо х,., , получим: 


Р(Г) 0 0 0 е!Т —е^Т 
г Ц е^Т .Р(Г) Р(Г) 0 0 р й е?"Т —е2^Т р 
ег" .Р(Г) е^.Р(т) РТ) 0 |’ "Пет ем "№ 
ез^Т .Р(Г) е?^Т .Р(Г) е^Т .Р(Г) Р(Г) е“"Т — е\^Т 


Матрица Е имеет размерность 8х4, а вектор Ё — размерность 8. 
Сформируем матрицу С и вектор №. Запишем ограничения (3) и (4) на управление и век- 
тор состояния в виде: 
и >И, 


Ик 2 Рак (10) 
Х к — Хы» 
Ха > т” 


для всех К = 0,3. 


Подставив выражения (9) для х,,, в систему неравенств (10), получим: 
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и >И, 
—ш >-Ц 
и >И, 
—м >-И 
и >И, 
—и, >-Ц 
и >И, 
и; > -И ть, 
е^' х, +Р(Т)и > Хи, 
—е^"х, - Р(Т)и, >-Х и, 
е^^Т х, +е^ТР(Г)и, +Р(Г)и, > Хи, 
—е”^" х, - е^"Р(Г)и, - Р(Г)и, > -Х ии, 
е^Т х, +е?^7Р(Г)и, + е^"Р(Г)и, + Р(Г)и, > Хи: 
—е°^Т х, — е?^Р(Т)и, -е^"Р(Т)и, - Р(Г)и, >-Х и, 
е“^Т х, +ез^"Р(Т)и, +е^^"Р(Г)и, + е^7Р(Т)и, + Р(Г)и, > Хи, 
—е“^Тх, - е—^"Р(Г)и, - е?^ТР(Г)и, - е^"Р(Ги, - Р(Г)и, >-Х д. 


тах / 


тах / 


тах / 


(11) 





Представляя неравенства (11) в матричном виде, получаем выражения для матрицы С 
размерности 16х4 и вектора А размерности 16: 














1 0 0 0 и 
1 0 0 0 не 
0 1 0 0 И 
0 о 0 0 и 
0 0 1 0 и 
0 0 Е 0 и 
0 0 0 1 й. 
д к 0 0 1 |. бы 
(Г) 0 0 0 Хи-е“"х 
т 0 оо Хо е”х, 
е”Р(Г) РТ) оо Хех, 
-етрт) РР) 0 о 0 Хы чех, 
е"Р(Т) е“РТ) РТ) 0 Хех 
-е"р(т) -е”рТ) РТ) 0 ине х, 
е”"Р(Т) е”“ТР(Т) е“Р(Т) РТ) х ах 
[-е?”р(т) -е?"Р(Т) -е^ТР(Т) -РТ). Хех, 


Результаты моделирования. Синтезируем систему управления электродвигателем постоянного 
тока ДПМ-25. На рис. 1 приведены графики переходного процесса в системе управления, эталон- 


ного переходного процесса и их производных, полученные при равенстве матриц 
К, =К, =К, = | , й ое | Выбор матрицы А на 4-м шаге со значительно большим 


первым диагональным элементом объясняется желанием обеспечить на последнем шаге мини- 
мально возможное отклонение положения от эталонного. 
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Рис. 1. Переходный процесс в системе управления электродвигателем 


Переходный процесс получен путём численного решения системы дифференциальных 
уравнений — (1), описывающих исследуемую систему с учётом электрической цепи якоря, с ша- 
гом равным 7/8. При синтезе учитывались ограничения на величину управляющего воздействия 


(напряжение на якоре) вида -—1 < ЦЕ) <1 и величину скорости вала двигателя вида -1 < «(Ё) <1. 


Для того чтобы гарантировать отсутствие перерегулирования вала двигателя, наложены ограни- 
чения с(Ё)>0. 


На рис. 2 изображён график ступенчатого управляющего воздействия. 
и(), В 


тах 





Рис. 2. Управляющее воздействие в синтезированной системе 


Заключение. На примере синтеза управления электродвигателем постоянного тока показано, 
что благодаря заданию желаемых свойств в виде эталонного переходного процесса и квадратич- 
ному критерию качества, характеризующего отклонение переходного процесса от эталонного, 
нахождение управления сводится к решению системы линейных уравнений с линейными ограни- 
чениями-неравенствами в реальном масштабе времени. 
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ОС МОТОК О1$СВЕТЕ СОМТВОЕ ТЕСНМТОЦЕ* 
А. К. Ау4туап, О. 1. Туе Кома, А. А. Мо!спапоу 


Тре (есптие оЁ пе @осгее сопго! оЁ а сопйпиоиз оБуЕСЕ [5 сопяегеа Бу {пе ехатр/е оРа ОС тоюг. Тре рго- 
рожей теёоа сап ятр/№у те ргосеадиге о! 5е/есипд те аезгей ргорегйе$ о! {те уппехгеа сопёго! зу${ет, апа 
гедисе {те сопЁго! Аепийсайоп Ю а Я&апаага ргоМет оп те [теаг едиайоп зузет зо/ивоп ит гедага ю {те 
ПедиаМу сопзатЁз аЁ еасй затр/е 5растд. Тйе аезгей ргорегиез оЁ {те зуппезгей зуйет аге зрестеа 
{ргоидВ {пе ге!егепсе апзепЕ. Тре диаагайс регогтапсе сшепоп #паЁ свагасейгез те аеиаНоп ог 1е {гап$!- 
оп ргосез$ Гот пе ге!егепсе аЁ те едиа!\ 5ерагаей тотепЁ оЁ ите [5 арр/еа. ТР арргоасв а!ои/5 а геа!- 
Ите сопёго! оЁ пе оБесЕ ий" ассоипЕ 0 Те гезтсНоп5 ог Бо Ю те сотго!тд асйоп, апа 1ю те рйазе сопго! 
зует соогтаЕе5. Трис, 1 15 ро55Ь/е юЮ ас!Геуе абжепсе о! оуег5ПооЕ апа озсИарИКу. 

Кеуигогаб: сопёго! асйоп, сопло! зучет, $ае иесюг, ОС тоюг. 





* Тпе гезеагсй {5 Чопе оп пете по. 01201259781 мййт {пе #ате ог {пе до\уегптепЕ {а$К ог {Ле Миху ог ЕдисаНоп апа $с1- 
епсе ог {Пе Кизчап Редеганоп т В&р. 
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